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We have measured the transverse thermal conductivity K xy 
of twinned and untwinned YBa2Cu307 single crystals as a 
function of angle 9 between the magnetic field applied par- 
allel to the CUO2 planes and the heat current direction, at 
different magnetic fields and temperatures. For both crystals 
we observed a clear twofold variation in the field-angle depen- 
dence of KxyiO) = —K xy (T, B) sin(2#). We have found that 
the oscillation amplitude K® y depends on temperature and 
magnetic field. Our results show that K® y = aBhi(l/(bB)) 
with the temperature- and sample-dependent parameters a 
and b. We discuss our results in terms of Andreev scattering 
of quasiparticles by vortices and a recently proposed theory 
based on the Doppler shift in the quasiparticle spectrum. 



To the main effects produced by the <i-wave symmetry 
of the superconducting order parameter one includes the 
influence of the extended quasiparticle (QP) states - as- 
sociated with a gapless structure (nodes) - on the mixed 
state of the superconductor in a magnetic field 0,^]. In 
a semiclassical way, the effect of a magnetic field can be 
taken into account by introducing a Doppler shift (DS) 
in the energy spectrum of the QP due to the superfluid 
flow around vortices. This idea was later developed to 
calculate the density of states, specific heat and thermal 
conductivity in the mixed state of d-wave superconduc- 
tors. Within this framework it is possible to explain 
the specific heat || and the thermal conductivity behav- 
ior in a magnetic field perpendicular to the CuC>2 planes 
at low enough temperatures 0. 

However, experiments on the field- angle dependence 
of the thermal conductivity with magnetic field applied 
parallel to the CuC>2 planes on YBa2Cus07 single crys- 
tals (Y123) were interpreted assuming Andreev scatter- 
ing (AS) of QP by vortices [js|-[i"cfl taking into account 
the d— wave symmetry of the order parameter. We think 
that this issue needs to be reconsidered experimentally to 
check whether there is clear evidence from thermal trans- 
port supporting either a DS effect on the energy spectrum 
of QP and/or AS. Furthermore, we note that there are 
no direct measurements published on the transverse ther- 
mal conductivity n xy for magnetic fields applied parallel 
to the CuC>2 planes, although the pioneer experiment by 
Yu et al. U provides an idea about the possible field- 
angle dependence of the transverse component. Nowa- 
days there is consent that K xy is a suitable property to 



check theoretical expressions of electronic thermal trans- 
port since it is free from non-electronic contributions. In 
this letter we report on the magnetic field B, field-angle 
6 and temperature T dependence of K xy for fields applied 
parallel to the CuC>2 planes. 

The magnetic field, field-angle and temperature depen- 
dence of K xy were studied on two Y123 single crystals. 
One sample was an optimally doped twinned single crys- 
tal with dimensions (length x width x thickness) 0.83 x 
0.6 x 0.045 mm 3 and critical temperature T c = 93.4 K 
p!l|,p2[ . The second sample was an untwinned Y123 sin- 
gle crystal with dimensions 2.02 x 0.68 x 0.14 mm 3 and 
T c = 88 K. The heat current J was along the longest 
axis of the crystal (in case of the untwinned (twinned) 
crystal it was parallel to its a- axis (a/6-axes)). The mag- 
netic field was applied parallel to the CuC>2 planes within 
a misalignment angle of ±0.5°. An in-situ rotation sys- 
tem enabled us the measurement of the conductivity as 
a function of the angle 9 defined between the applied 
field and the heat current direction. The used tempera- 
ture gradient V ' X T < 300 K/m was kept constant along 
the heat current direction x during all scans in order to 
calculate the transverse component of the thermal con- 
ductivity using the relation 



K xy (T,B,9) = k xx (T,B, 
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In this equation we used the experimental fact that the 
oscillation amplitude of the longitudinal component of 
the thermal conductivity k xx - which has been shown 
to have fourfold symmetry at T/T c < 0.26 is 
much smaller than the longitudinal temperature gradi- 
ent | V X T| |l5[ ). Therefore n xy shows the same field-angle 
dependence as V y T. 

The longitudinal and transverse temperature gradients 
were measured using a previously field- and temperature- 
calibrated type E thermocouples with a dc-picovoltmeter 
Jhlj . We obtained similar results in both crystals for the 
field-angle and field dependence of K xy . The tempera- 
ture dependence of K xy as well as its absolute value de- 
pends, however, on the sample. This result appears to be 
related to different impurity scattering rates. The mea- 
sured behavior of the longitudinal thermal conductivity 
as a function of temperature indicates a stronger impu- 
rity scattering for the untwinned crystal probably related 
to its oxygen deficiency. 
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FIG. 1. (a) Angle dependence of the transverse temper- 
ature gradient at a fixed field of 8 T and at different tem- 
peratures below T c for the twinned crystal. Continuous lines 
are fits to the function sin(2(9). (b) Angle dependence of the 
transverse conductivity at T — 29.8 K at different applied 
fields parallel to the Cu02 planes for the twinned crystal. If 
the heat current direction is along +x and a positive 90° an- 
gle is along +y, then there is a positive transverse thermal 
gradient for B at +45°. 

Figure |](a) shows the transverse temperature gradient 
V y T as a function of the angle 9 at different temperatures 
at a constant field B = 8 T. A clear twofold pattern is 
measured as the magnetic field is rotated parallel to the 
CuC>2 planes. Due to the enhancement of the inelastic 
scattering rate of QP with increasing temperature the 
angle dependence cannot be resolved above 70 K within 
the experimental error. The background in each curve 
of Fig. ||(a) is due to the experimental misalignment in 
the position of the sensors of the thermocouple used to 
measure V y T since its temperature dependence follows 
roughly that of \7 X T (the transverse thermocouple mea- 
sures partially the longitudinal gradient on the sample 
as well). Taking into account the measured longitudi- 
nal conductivity k xx at B = we estimate a maximum 
misalignment of 80 fmx along the heat current direction 
between the ends of the transverse thermocouple sen- 
sors. Thus, since the length and width of both crystals 



are larger than the misalignment of the transverse ther- 
mocouple, we can consider our system as a probe of the 
true transverse signal, in contrast to Ref. ||. 

A proof of this assumption is that within our experi- 
mental resolution the curves keep their background when 
the magnetic field is changed. This result allows us to 
subtract it from the curves and to consider only the os- 
cillation amplitude K® y (T, B) which contains the physical 
information of the transverse response of the system. The 
angular dependence of K xy obtained in this way is shown 
in Fig. 0(b) for different values of the magnetic field. We 
think that this is the first time that such a symmetry is 
unambiguously established for n xy ; it follows a field- angle 
dependence proportional to sin(20). Based on the pub- 
lished work we assume that the twofold field-angle depen- 
dence reflects the symmetry of a d— wave order parame- 
ter. Note that the amplitude of the angle dependence on 
the transverse gradient is a relatively larger effect than 
in the longitudinal case (e.g. for the twinned crystal at 
B = 8T and T — 13.8 K, V y T = 2.1 K/m for a longitu- 
dinal gradient of V X T = 241 K/m; the amplitude of the 
oscillation of the longitudinal component remains of the 
order of one percent or less of V X T). This is the reason 
for the results obtained in the experiment of Yu et al. 
H where only a twofold pattern was measured, not the 
fourfold one expected for the field-angle dependence of 
the longitudinal thermal conductivity |p|,fL0| . A compari- 
son between these two effects will be reported elsewhere 
II- 

The observed angular dependence of the thermal con- 
ductivity in Refs. p|-|i"o[ was interpreted in terms of AS 
of QP by vortices screening currents using the 2D ex- 
pression of the Bardeen-Richayzen-Tewordt model and 
assuming a d— wave order parameter (l^^] : 

(2) 

where, in general, T(B 7 p, T) is the relaxation rate given 
by the sum of scattering rates of QP by impurities 
l/ro(p), by phonons l/r P (p, T), by QP 1/t(B, p, T) and 
AS by vortex currents 1/t v (B, p, T) a and j3 de- 
note the x or y directions on the plane of the sample 
and v ga is the group velocity along the a direction. All 
other symbols have the conventional meaning and are 
described in Ref. Q in detail. With similar parameters 
as used in f|,[l(J, e.g. B%% = 500 T, Ginzburg-Landau 
parameter k = 100 and a zero-temperature energy gap 
Ao(0) = 20 meV we obtain a sin(20) angle dependence 
for K xy as the experiments show. In what follows we com- 
pare our experimental results and those from literature - 
for magnetic fields applied parallel to the Cu02 planes - 
and discuss to which extent the available data can help to 
differentiate between the AS and DS contributions pro- 
posed by the existing models. 
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FIG. 2. Temperature dependence of the normalized oscil- 
lation amplitude ft!L at different applied fields for the un- 
twinned crystal. The two continuous curves are obtained with 
(2) with the parameters mentioned in the text and normalized 
by k° h (55K,8.5T). K° xy (U K,8.5 T) = 0.084 W/Km. 

(1) Angular dependence of n xy (0): We note that the 
observed angle dependence of V y T agrees with the ex- 
pected sign if AS affects the thermal transport. For ex- 
ample, if B is at +45° the measured V y T, see Fig. [|, 
indicates a higher T at +90°, i.e. there is an excess of 
QP flowing along the field relative to those moving in the 
—45° nodal direction. 

Recent studies pointed out that a DS in the excitation 
spectrum of QP would produce also a twofold pattern in 
the field- angle dependence of k xv respectively [lfjJlT]]. In 
particular, Won and Maki [ p^Tq ] obtained the following 
expression for the transverse conductivity: 
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that applies for T -C T c , B -C B C 2 and in the clean 
limit To S> h/A, where A(r) ~ Ao, this last an impurity 
independent order parameter. Here T = t q ~ 1 /2 repre- 
sents the impurity scattering rate, Ko = hkgTn/3Am, 
v = vp, v' = VFX a /X c , n the density of QP and m their 
mass, b = hvv'e/AA 2 . The clean-limit expression (|3|) 
gives a simple field dependence for K xy with both, a and b, 
temperature and sample dependent parameters in which 
the scattering rate T of the QP is included. Within a 
simple picture the sign of the measured V y T seems to 
disagree with that expected from the DS: an excess of 
QP should be in the nodal direction perpendicular to the 
applied field. However, the DS affects both the carrier 
density as well as the scattering rate. This last dominates 
K xy at high enough temperatures |l8| , thus the QP move 
easier parallel to the magnetic field and the sign in (||) 
turns out to be the correct one. The temperature depen- 
dence of K xy in the clean limit is still under development. 



Although preliminary results [18| indicate that the main 
arguments in Eq. (0) may remain, we should take @ as 
a trial function to compare with our results. 
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FIG. 3. Magnetic field dependence of the normalized os- 
cillation amplitude at different temperatures for the twinned 
(□) and untwinned (•) crystals. The dashed and dotted lines 
are fits of the data to the function B ln(l/(feB)) with 6 as a 
free parameter. The continuous line is obtained from Eq. (^) 
assuming Andreev and impurity scattering with the same pa- 
rameters as in Ref. |]]. 

(2) Temperature dependence of n xy : This depen- 
dence is presented in Fig. || for different magnetic fields. 
We note that every point in this figure represents the 
measurement of a complete curve as shown in Fig. |l| 
and that the experiments below 20 K have to be per- 
formed in a field-cooled procedure to avoid pinning effects 
]l0|,[ll]]. The results show a maximum in K° xy at a field- 
independent temperature of T ~ 35 K. n xx (T) shows a 
maximum at a similar temperature, also for fields applied 
parallel to the c-axis of the crystal. These facts indicate 
a common origin for the observed maxima: a competi- 
tion takes place between a decreasing inelastic scattering 
rate due to the opening of the energy gap and the de- 
crease in the number of QP with decreasing temperature 
]l9| , p0[ . Therefore, to fit the observed T— dependence we 
need the rates l/rp(p, T) and l/r(B,p,T) to include in 
r, see Eq. (||). As seen in Fig. 2, with the above men- 
tioned parameters the AS model and without l/rp(p, T) 
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and 1/t(B, p, T) does not provide a satisfactory fit of the 
data. Unfortunately, the exact expressions as well as the 
values of the necessary parameters of the rates are not 
known with certainty to allow a confidential result from 
the fits. Then and since Eq. (Q) is strictly valid at low 
T, from the observed T— dependence it is not straight- 
forward to conclude which of the proposed mechanisms 
predominates. 
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FIG. 4. Temperature dependence of the parameters a(«) 
and &(o) obtained from the fits to Eq. (^) (see Fig. ^ for the 
untwinned (a) and twinned (b) crystals. The parameter a is 
obtained taking the absolute value of k® v (T). 

(3) Field dependence of n xy : This dependence is pre- 
sented in Fig. U for the two samples at different constant 
temperatures. Equation (||) describes satisfactorily the 
field dependence of K xy from the lowest measured temper- 
ature T ~ 13 K. The temperature dependence of the fit 
parameters a and b are presented in Fig. [|. We note that 
a and b show a similar T— dependence for the two sam- 
ples but different absolute values. From the parameter a 
obtained at T < 20 K we calculate the scattering time 
to ~ 0.34(2.9) ps for the untwinned (twinned) crystal as- 
suming A = A o (0) = 20 meV and v ~ 5v' ~ 2 x 10 5 m/s. 
With these values we obtain a reasonable agreement for 
the b parameter of the untwinned crystal. Note, how- 
ever, that the predicted scattering-rate dependence in b 
is not only too small in comparison to the factor 10 differ- 



ence between samples, but is qualitatively different, i.e. 
b increases with V f2lf| . Assuming Andreev and impurity 
scattering rates in (0) we cannot reproduce satisfacto- 
rily the observed field dependence of Fig. || in a broad 
temperature range even varying generously the fitting 
parameters. In Fig. || we show the results of this cal- 
culation using similar parameters as indicated above. At 
high T both models provide similar good fits to the data. 

(4) Field and angular dependence of k xx : In the tem- 
perature region of our work k xx decreases with field (see 
for example Fig. 15 in |11|). This effect as well as the an- 
gular fourfold pattern [ 10 1 appear to be accounted for by 
the AS mechanism and at odds with the DS predictions 
since in the clean limit (dirty limit in [||) a DS increases 
K x , r and gives the incorrect sign for the fourfold pattern 
fll6[ . However, two new approaches can provide an expla- 
nation of the results within the DS effect. First, high T 
and the influence of scattering processes revert the field 
dependence of k xx (as shown already for the superclean 
limit in ||) and the sign of the fourfold pattern [Q. Sec- 
ond, we note that experiments measure the total conduc- 
tivity in which phonons also contribute as heat carriers, 



i.e. k xx {T, B) = K e J x (T,B) 



. P h 



(T, B) . In general it is 



assumed that kP x \ is field independent. This is not nec- 
essarily true. Note that none of the separation methods 
used in literature is accurate enough to prove that kP x \ 
is strictly field independent. Due to the phonon-clcctron 
interaction a decrease of kP xx of a few percent, for fields 
of a few teslas and parallel to the planes, is possible. In 
this case the DS may decrease the total conductivity with 
field due to its influence to the phonon contribution |Q . 
This contribution would provide also the correct sign of 
the fourfold pattern. 

(5) Angular dependence of the specific heat: Recent 
specific heat results indicate the absence of a fourfold 
symmetry in its angular dependence, expected if the DS 
changes the QP density of states [£2| . It is still unclear, 
however, how large this effect should be Jl6| as well as the 
influence of extra effects that may not allow an accurate 
test of the expected fourfold symmetry in the specific 
heat §§. 

In summary, a clear twofold angle dependence has been 
obtained in the whole measured temperature and field 
range for the transverse thermal conductivity. Its oscil- 
lation amplitude n xy follows a(T)B ln(l /b(T)B) with the 
sample-dependent parameters a and b. Rigorously speak- 
ing none of the available models explain satisfactorily all 
the details of the experimental data. The expected influ- 
ence of the DS casts doubts whether the AS mechanism 
is the only and appropiate one to understand K xy (T, B). 
We think that a further development of both models is 
necessary before a clear answer can be given. 
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